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Optical  time  domain  ref lectometry  (OTDR)  has  been  suggested  as  a  tool 
for  detecting  intrusions  in  a  fiber.  This  is  based  on  the  knowledge 
that  the  OTL)R  is  very  sensitive  in  locating  and  measuring  fiber  per¬ 
turbations,  including  loops,  bubbles,  and  sinusoidal  perturbations  of 
optic  axis  and  core  diameter. 

The  program  goal  is  to  develop  the  OTDR  for  intrusion  monitoring  and  to 
establish  its  ultimate  performance  limits.  An  additional  objective  is 
to  identify  and  catalog  OTDR  signatures  of  different  potential 
intrusions . 


'Theoretical  considerations  indicated  that  maximum  perturbation  sensi¬ 
tivity  of  the  OTDR  is  achieved  by  increasing  the  source  power  until  the 
threshold  for  stimulated  effects  is  reached.  This  maximizes  sensi¬ 
tivity  regardless  of  the  relative  shapes  of  the  probe  pulse  and  the  per¬ 
turbation  response.  Tailoring  the  probe  pulse  shape  to  be  a  matched 
filter  for  the  perturbation  results  in  a  further  increase  in  sensi¬ 
tivity  at  the  cost  of  greatly  increased  complexity  in  the  source. 

!  Several  characteristics  of  perturbation  signatures  have  been  found.  Per¬ 
turbations  which  are  much  smaller  in  spatial  extent  than  the  probe  pulse 
are  characterized  by  the  shape  of  the  probe  pulse.  If  the  probe  pulse 
width  is  much  smaller  than  the  spatial  extent  of  the  perturbation,  the 
OTDR  response  is  the  impulse  response  of  the  perturbation. 

Experiments  on  splices  indicate  that  it  is  very  difficult  to  splice 
together  two  waveguides  so  perfectly  that  the  OTDR  is  unable  to  detect 
the  splice.  In  every  observed  case,  an  insertion  loss  was  measured  and 
conversion  of  bound  mode  power  to  leaky  modes  was  observed. 

Experiments  on  measuring  the  excess  loss  of  a  single  fiber  loop  indi¬ 
cated  that  the  OTDR  technique  is  capable  of  detecting  loops  as  large  as 
8-cm  in  diameter.  Larger  diameter  loops  can  be  located  if  the  system 
signal -to-noise  ratio  is  improved. 

Large  as  well  as  small  amplitude  microbends  in  a  step-index  waveguide  can 
easily  be  located  and  measured  with  the  OTDR. 

A  singj.e_microbend  less  than  10  pm  in  amplitude  and  a  periodic  micro¬ 
bend  ^A03periods)  of  Amplitude  less  than  2  pm  were  detectable  with  a 

computerized  OTDR  test  set,  '  \ 
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SECTION  1 


INTRODUCTION  AND  SUMMARY 


Wideband  communications  systems  of  the  near  future  will  use  optical  fibers 
as  the  transmission  medium  instead  of  conventional  conducting  metal  waveguides. 
The  limiting  factors  of  conventional  systems  are  high  attenuation,  high  costs, 
and  large  size  and  weight.  By  contrast,  the  wideband  frequency  capability  of 
optical  fibers,  the  development  of  very  low  loss  optical  fibers  (<1  dB/km),  the 
rapid  development  of  associated  optical  circuit  elements  (such  as  connectors, 
optical  sources,  detectors,  optical  couplers,  switches,  and  multiplexers),  and 
the  promise  of  eventual  low  costs  for  these  key  elements  are  expected  to  result 
in  widespread  use  of  optical  fiber  communication  systems.  However,  hostile 
intrusion  into  these  optical  communication  links,  particularly  into  DOD  links, 
will  require  immediate  attention  and  studies,  along  with  a  method  for  detecting 
these  intrusions,  before  hard  standards  for  secure  fiber  communication  links  can 
be  set. 

Recently,  a  new  concept  for  detecting  intrusions  using  an  optical  time 
domain  ref lectometer  (OTDR)  has  been  suggested.  This  suggestion  is  based  on  the 
knowledge  that  the  OTDR  is  a  sensitive  tool  for  locating  and  measuring  fiber 
perturbations  including  loops,  bubbles,  and  sinusoidal  perturbations  of  optic 
axis  and  core  diameter. 

The  great  advantage  of  OTDR  monitoring  is  that  the  intrusion  point  can  be 
quickly  determined  and  corrective  action  applied  at  that  point.  Although  the 
technique  currently  employed  for  monitoring  intrusion  —  by  monitoring  the  change 
in  transmitted  light  intensity  or  modal  characteristics  of  the  light  transmitted 
from  station  to  station  —  is  adequate  for  detecting  an  intrusion,  it  does  not 
appear  capable  of  providing  information  that  can  be  used  to  localize  the  point 
of  intrusion.  Thus,  OTDR  monitoring  provides  an  attractive  complement  to  trans¬ 
mission  monitoring. 

There  are  various  methods  for  tapping  power  from  an  optical  fiber  communi¬ 
cation  system  (as  there  are  from  copper  communication  links).  Periodic  deforma¬ 
tion  of  the  fiber,  microbending  of  fibers,  proximity  coupling,  etc.  are  methods 
suitable  for  tapping  power  from  an  optical  fiber.  These  intrusion  methods 
result  in  power  loss  from  the  fiber  and  are  detectable  with  the  OTDR.  In  addi¬ 
tion  to  localizing  the  position  of  intrusion,  it  is  desirable  to  characterize 
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these  intrusions  with  "signatures”  using  the  OTDR.  It  may  be  possible  to  relate 
different  types  of  intrusions  to  different  signatures.  It  is  one  of  the  goals 
of  this  research  program  to  determine  whether  a  high-sensitivity  OTDR  can 
resolve  signatures  characteristic  of  the  intrusion. 

An  additional  potential  advantage  of  determining  intrusion  signatures  is  to 
compare  these  signatures  with  those  occurring  from  environmental  events  affect¬ 
ing  the  fiber.  Environmental  events  include  random  or  steady  winds,  rapid  tem¬ 
perature  changes  resulting  in  the  fiber  rubbing  in  its  cabling,  strain-induced 
loss  in  the  fiber,  and  microbending  losses.  The  resulting  fiber  losses  may 
cause  frequent  false  intrusion  alarms  if  conventional  monitoring  is  used  alone 
and  could  result  in  substantial  link  down-time.  If  the  signatures  of  environ¬ 
mental  events  are  found  to  be  considerably  different  from  those  of  real  intru¬ 
sions,  the  down-time  of  the  communication  link  can  be  decreased  significantly  by 
shutting  down  the  link  only  when  an  intrusion  signature  appears. 

The  purpose  of  this  program  is  to  develop  the  OTDR  for  intrusion  monitoring 
and  to  establish  its  ultimate  performance  limits.  A  primary  objective  of  this 
program  is  to  determine  both  the  ultimate  resolution  and  the  ultimate  sensitiv¬ 
ity  of  the  OTDR.  An  additional  objective  of  this  effort  is  to  identify  and 
catalog  OTDR  "signatures"  of  different  potential  intrusions.  Since  each  type  of 
potential  intrusion  affects  the  bound  modes  of  the  fiber  in  a  different  way,  the 
OTDR  characterization  of  each  intrusion  is  unique.  Consequently,  a  catalog  of 
intrusion  signatures  can  be  made. 

Major  accomplishments  of  this  program  are  listed  below: 

•  Theoretical  considerations  of  the  maximum  sensitivity 
of  the  OTDR  technique  to  the  presence  of  perturbations. 

o  Modifications  and  improvements  of  the  OTDR  test  set  to 
facilitate  measurements  of  perturbations. 

«  Detection  of  the  presence  of  a  single  loop  in  an  optical 
fiber. 

•  Detection  of  the  presence  of  fiber  splices  and  raicrobend 
transducers . 
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Measurement  of  the  insertion  loss  of  a  single  raicrobeand. 

Measurements  of  the  diameter  variations  of  a  fiber  and  its 
correlation  to  the  OTDR  signature. 

Investigation  of  the  effects  of  a  periodic  microbend  trans 
ducer  on  fiber  loss. 


SECTION  2 


REVIEW  OF  THE  OTDR  TECHNIQUE 

The  method  of  measuring  waveguide  attenuation  with  the  optical  time  domain 
ref lectometer  (OTDR)  is  often  referred  to  as  the  backscattering  method.  It 
measures  the  time  rate  of  change  of  light  Rayleigh  backscattered  from  an  optical 
pulse  propagating  along  a  waveguide.  The  rate  of  change  is  'tted  to  an  expo¬ 
nential,  exp(-2ax),  where 


and  c/n  is  the  group  velocity.  Attenuation  A  is  defined  by 


A  =  4.34  ax 


If  the  distance,  x,  is  measured  in  kilometers,  attenuation  has  units  of  dB/km. 

This  method  has  many  advantages.  In  particular,  it  is  unique  in  providing" 
information  on  the  length  dependence  of  the  loss.  It  is  non-destructive,  accur¬ 
ate,  and  highly  sensitive  to  small  perturbations.  A  singular  advantage  of  the 
technique  is  that  It  normally  requires  access  to  only  one  end  of  the  waveguide 
in  order  to  perform  loss  measurements.  This  last  feature  obviously  makes  it 
extremely  attractive  for  use  in  field  Installations. 

However,  the  method  has  two  disadvantages:  the  exact  modal  distribution  of 
the  backscattered  power  is  unknown,  and  an  ambiguity  can  result  when  particu¬ 
larly  nonuniform  waveguides  are  measured.  The  ambiguity  can  be  resolved  by 
measuring  the  waveguide  from  both  ends  or  by  operator  interpretation  of  the 
data. 

*  The  backscattering  method  permits  the  measurement  of  nearly  every  waveguide 
parameter.1  Included  in  the  measurement  capabilities  are:  total  attenuation 
coefficient,  scattering  coefficient,  absorption  coefficient,  differential  mode 
attenuation,  waveguide  length,  bandwidth,  location  of  defects,  faults,  connec¬ 
tors,  splices,  and  the  effect  of  the  external  environment  on  waveguide  proper¬ 
ties  (including  introduction  of  lossy  coatings,  strain,  and  ionizing 
radiat ion) . 


11 


The  general  experimental  arrangement  used  in  the  OTDR  is  shown  in  Figure  1. 
Qualitatively,  light  from  a  pulsed  source  is  coupled  into  a  waveguide  using  an 
optical  directional  coupler.  As  the  pulse  propagates,  it  encounters  a  string  of 
absorption  and  scattering  sites  which  attenuate  the  beam  (see  Figures  2(a)  and 
(b);  each  scattering  site*  affects  the  forward-traveling  pulse  in  a  unique  way. 
However,  each  site  converts  some  of  the  forward  traveling  pulses  into  a  reverse 
traveling  pulse  (Figure  2(c)).  That  portion  of  the  reverse-traveling  pulse 
which  is  trapped  in  the  guided  modes  will  be  directed  back  to  the  directional 
coupler.  As  the  backscattered  pulse  propagates,  it  is  again  attenuated  by  the 
scattering  and  attenuation  sites  it  encounters.  When  the  resultant  signal 
(Figure  2(d))  reaches  the  directional  coupler  (Figure  1)  it  is  routed  to  a 
photodetector.  The  detected  signal  is  amplified  and  then  electronically 
processed  before  being  displayed  and  recorded. 

A.  ANALYSIS  OF  AN  UNPERTURBED  WAVEGUIDE 

The  backscattering  method  can  be  analyzed  quantitatively  by  considering  a 
pulse  of  light  propagating  along  a  waveguide.  The  forward-  raveling  pulse  suf¬ 
fers  a  position-dependent  total  attenuation  determined  by  the  total  attenuation 
coefficient,  Yft>  which,  in  general,  depends  upon  position  and  mode  group 
excitation.  This  coefficient  consists  of  two  terms,  a  scattering  term, 

Yfs,  and  an  absorption  term,  Yfa* 

The  scattering  terra,  Yfg,  in  low-loss  w;weguldes  is  due  predominantly 
to  Rayleigh  scattering  from  two  main  sources:  density  fluctuations  and  concen¬ 
tration  fluctuations.  Density  fluctuations  are  thermal  fluctuations  of  the 
glass  atoms  which  are  frozen  in  during  waveguide  fabrication.  The  scattering 
coefficient,  Y^d,  due  to  density  fluctuations,  is  related  to  the  iso¬ 
thermal  compressibility,  8,  through  the  expression, 


*Scattering  sites  are  defined  here  to  mean  any  site  which  converts  a  portion  of 
the  f orward-t ravel Lng  wave  into  a  reverse-traveling  wave.  According  to  this 
definition,  scattering  sites  include  dielectric  discontinuities  resulting  in 
Fresnel  reflection,  Rayleigh  scattering  sites,  and  point  defects  that  cause 
reflect  ions. 
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WAVEGUIDE 
UNDER  TEST 


Figure  1. 

A  block  diagram  of  the  OTDR.  A  pulse  of  light 
from  the  source  is  coupled  into  the  waveguide 
under  test  by  an  optical  directional  coupler. 
The  directional  coupler  also  routes  the  back- 
scattered  signal  to  a  photodetector.  The 
detected  signal  is  amplified  and  electronically 
processed  before  recording  and  display. 
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1 1567  -3  i. 


Figure  2. 

The  backscatteriag  process  la  a  waveguide  of  length  L.  The 
waveguide  contains  Innumerable,  randomly-distributed,  micro¬ 
scopic  scattering  and  absorption  sites,  (a).  As  the  probe 
pulse  propagates  (b),  it  is  attenuated  by  these  sites.  At 
point  zc  the  pulse  is  backseat tered.  The  backscattered  signal 
is  again  attenuated  by  the  scattering  and  absorption  sites,  (c). 
The  resultant  waveform  is  seen  in  (d).  It  has  been  plotted  on  a 
new  abscissa,  z’,  to  show  that  this  waveform  attenuates  at  twice 
the  rate  of  the  probe  pulse  if  we  plot  it  against  the  length  of 
the  waveguide. 


(3) 


\d '  ^  ("2  “ 11  kT6 

where  n  is  the  refractive  index,  k  is  Boltzmann's  constant,  T  is  the  transition 
temperature  at  which  the  density  fluctuations  are  frozen  in,  and  X  is  the  wave¬ 
length.  For  fused  silica,2  T  is  ~800°K  and  8  is  ~7  x  10” 11  m3/J. 

If  a  glass  waveguide  is  fabricated  by  introducing  an  oxide  other  than  Si02 
to  form  a  region  of  slightly  different  index,  concentration  fluctuations  of  the 
dopant  oxide  will  also  cause  scattering.  The  scattering  coefficient,  Ysc, 
due  to  concentration  fluctuation,  is 


where  (Ac)2  is  the  mean-square  concentration  fluctuation,  and  V  is  the  volume 
over  which  the  concentration  fluctuation  occurs.  For  high  silica  content  glass 
waveguides  this  scattering  mechanism  accounts  for  about  25%  of  the  total  scat¬ 
tering.2  Neglecting  any  other  scattering  mechanisms,  the  total  scattering 
coefficient,  Yfs,  is 


TfS  =  vSd  +  V 


(5) 


The  waveguide  is  characterized  by  this  coefficient  and  Yfa»  the  absorption 
coefficient: 


ft 


Yfs  + 


fa 


(6) 


where  differential  modal  attenuation  is  neglected. 

The  experiment  shows  that  as  a  pulse  propagates  down  a  waveguide  the  change 
in  power,  AP,  in  an  incremental  distance,  Ax,  is  proportional  to  the  power  at 
that  point.  Thus, 
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APft(x)  =  -yft(x)Pft(x)Ax  . 


(7) 


The  minus  sign  indicates  that  the  power  decays  as  the  pulse  propagates.  In  the 
limit  as  Ax  ♦  0,  Equation  (7)  can  be  rewritten  as 


Pft(z)  =  P0exp 


Z 

J  Yft(x)dx  . 
o 


(8) 


The  backseat tered  power,  Pbs,  is  some  fraction,  f(z),  of  the  power  at 
position  z: 


Pbs(z)  =  f(z)Pft(z)  .  (9) 

As  the  backscattered  signal  propagates  toward  the  source,  it  suffers  a  position- 
dependent  total  attenuation,  Ybt(z)*  Consequently,  the  power  backscat¬ 
tered  from  the  pulse  at  position  z  that  eventually  finds  its  way  back  to  the 
detector  is 


P(z)  =  p0f(z)  exP 


z 

( |T{tU)  +  v<x) 

*0 


> 


(10) 


where  P0  Is  the  initial  peak  pulse  power. 

The  better  waveguides  made  today  exhibit  practically  no  attenuation  other 
than  Rayleigh  scattering  from  microscopic  density  fluctuations  frozen  in  during 
fabrication.3  For  such  a  waveguide,  the  backscatterlng  function  is  a  constant,4 


f(z) 


I  XNAll  .  _J2_ 

8  Ts  2  a  +  1 
n„ 


(11) 
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Here  NA  is  the  numerical  aperture,  nQ  is  the  axial  refractive  index,  yg  is  the 
scattering  coefficient  (where  we  assume  Yfs  *  Ybs  =  Ys)>  an(1  a  is 
the  refractive  index  profile.  This  is  defined  in 


2 

n2(r)  - 

n0 

2 

n2 

1  -  2A  ( 

if 

r  <  a 
r  >  a  , 


(12) 


where  r  is  the  radial  coordinate,  a  is  the  core  radius,  n2  is  the  refractive 

index  of  the  cladding,  and  A  =*  (n2  -  n2)/2n2  is  related  to  the  index  difference 

o  2  o 

between  core  and  clad.  If  the  waveguide  attenuation  is  the  same  for  pulses 
traveling  in  both  directions,  the  attenuation  coefficients  are  identical;  thus 


Tft(x)  -  Ybt(*>  5  Yt 


(13) 


and 


y 


t 


y_  +  y. 


In  this  case,  Equation  (10)  reduces  to  Beer's  law: 


(14) 


P(z) 


,  y  (NA)V, 

3.  a  s  0 

8  a  +  1  2 


z 


(15) 


The  OTDR  signature  of  such  a  commercially  available  waveguide  is  shown  in 
Figure  3.  The  word  signature,  used  to  describe  the  OTDR  characterization  of  a 
waveguide,  emphasizes  that  each  waveguide  is  unique  and  exhibits  a  unique 
characterization. 

Three  features  are  apparent  In  the  signature  of  Figure  3.  First  is  the 
Fresnel  reflection  which  occurs  at  the  air-glass  interface  as  the  probe  pulse 
enters  the  waveguide.  This  reflection  is  generally  undesirable.  Next  is  a 
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RELATIVE  BACKSCATTERED  POWER 
ARBITRARY  UNITS 


POSITION  ALONG  WAVEGUIDE,  m 

Figure  3. 

,  An  OTDR  characterization  of  a  low-loss  commercially 

*  available  waveguide.  The  measurement  was  made  at  0.9  ym. 

Note  the  Fresnel  reflections  which  occur  at  each  end  of 
the  waveguide  and  the  exponential  decay  of  the  Rayleigh 
scattered  light  between  them. 

i 
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region  of  exponential  decay  of  the  Rayleigh  scattered  light.  Finally,  we  see 
the  Fresnel  reflection  which  occurs  as  the  pulse  leaves  the  far  end  of  the  wave¬ 
guide.  The  amplitude  of  the  reflection  from  the  far  end  can  be  used  to  deter¬ 
mine  the  quality  of  the  cleave,  as  well  as  the  angle  that  the  cleaved  surface 
makes  with  the  optic  axis.  The  length  of  the  waveguide  can  be  determined  by 
measuring  the  time  delay  between  the  two  Fresnel  reflections  and  multiplying  by 
half  the  group  velocity  of  light  in  the  core. 

B.  ANALYSIS  OF  A  PERTURBED  WAVEGUIDE 


The  previous  discussion  assumed  that  the  waveguide  was  perfect,  but  most 
waveguides  are  imperfect.  To  determine  the  OTDR  response  of  an  imperfect  wave¬ 
guide,  the  probe  pulse  is  labeled  as  Ap(ct-z)  and  the  unperturbed  waveguide 
response  as  b(z),  where  A  is  a  constant.  Assuming  the  waveguide  response 
remains  linear  when  a  small  disturbance  d(z)  is  added,  the  backscatter  power 
B(t)  can  be  written  as 


B(t)  -  Aj[b(z)  +  d(z)]  p(ct-z)  dz  .  (16) 

Equation  (16)  is  just  the  sum  of  convolutions, 

B  =  Ab*P  +  Ad*P  ,  (17) 

where  the  symbol  (*)  denotes  the  convolution  operation.  We  see  from 
Equation  (17)  that  the  OTDR  signature  consists  of  two  terras:  the  unperturbed 
response  and  the  response  due  to  the  disturbance. 

There  are  three  classes  of  perturbations  of  special  interest:  a  very  nar¬ 
row  probe  pulse,  a  point  defect,  and  defects  where  the  spatial  extent  of  the 
probe  pulse  is  approximately  equal  to  that  of  the  defect.  Since  the  operation 
of  convolution  is  commutative,  the  first  two  cases  are  formally  equivalent.  In 
the  first  case,  the  probe  pulse  only  samples  the  perturbation  over  an  incre¬ 
mental  distance  much  less  than  the  width  of  the  perturbation.  In  this  case,  we 
have 
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(18) 


f*p  *  A  J d(z)  p(ct-z)  dz 
=  A  f  d(z)  <5  (ct-z)  dz 
=*  A  d(ct) 

Equation  (18)  states  that  the  backscattered  power  is  a  faithful  reproduction  of 
the  perturbation  response  (scaled  by  the  amplitude  of  the  probe  pulse)  when  the 
perturbation  is  much  larger  than  the  width  of  the  probe  pulse. 

In  the  second  case,  the  width  of  the  probe  pulse  is  much  greater  than  the 
perturbation,  and  the  convolution,  d*p,  becomes 

d*p  =  A  p(ct)  .  (19) 

Equation  (19)  states  that  the  backscattered  power  is  a  faithful  reproduction  of 
the  probe  pulse.  Hence,  a  narrow  probe  pulse  yields  the  impulse  response  of  the 
perturbation.  On  the  other  hand,  a  point  defect  can  be  identified  by  its 
resemblance  to  the  probe  pulse. 

In  the  last  case,  where  the  perturbation  is  approximately  the  same  width  as 
the  probe  pulse,  details  of  their  shapes  become  important.  In  general,  the 
response  will  not  be  a  faithful  reproduction  of  either  the  probe  pulse  or  the 
perturbation,  but  will  be  wider  than  either. 

This  analysis  indicates  the  tremendous  potential  of  this  technique  for 
studying  defects  in  waveguides.  Point  defects  can  be  located  with  the  OTDR  and 
then  examined  to  determine  the  source  of  the  defect.  Extended  defects  affect 
bound  mode  power  in  different  ways;  therefore,  examining  the  OTDR  signature 
provides  clues  about  the  origin  of  the  defect. 


20 


SECTION  3 


THE  OTDR  EXPERIMENTAL  SETUP 


Figure  4  shows  the  block  diagram  of  an  automated  OTDR  used  in  this  program. 
One  of  the  main  features  of  this  setup  is  the  use  of  a  polarization-isolation 
directional  coupler  based  on  a  Glan-prism  polarizer.  This  coupler  is  used  in 
the  OTDR  since  it  reduces  mount  reflection,  permits  launch  mode  control,  and  is 
a  non-destructive  power  launching  method.  Light  from  an  injection  laser  is 
coupled  into  a  fiber  pigtail.  (The  pigtail  facilitates  source  changes  and  pro¬ 
vides  a  more  nearly  circularly  symmetric  launched  spot.)  The  light  output  from 
the  pigtail  is  collimated  by  the  microscope  objective,  L, .  The  collimated  beam 
should  be  just  large  enough  to  completely  fill  the  aperture  of  the  lem,  L-,. 

The  Glan-prism  polarizer  transmits  one  polarization  and  dumps  the  power  in  the 
other  polarization.  For  discussion  purposes,  we  assume  that  the  horizontal 
polarization  is  transmitted  and  the  vertical  is  lost.  A  second  microscope 
objective,  L2>  focuses  light  onto  the  waveguide  under  test.  (Maximum  coupling 
efficiency  is  achieved  if  this  lens  has  a  numerical  aperture  identical  to  that 
of  the  waveguide.)  The  Fresnel  reflection  from  the  cleaved  input  endface  of  the 
waveguide  under  test  remains  horizontally  polarized.  As  a  result,  this 
reflected  light  is  not  directed  toward  the  detector;  instead  it  passes  through 
the  prism  polarizer  and  is  lost.  The  backscattered  power  is  collimated  by 
lens  L2«  Vertically  polarized  light  is  directed  toward  the  detector  by  the 
polarizer.  (Horizontally  polarized  light  is  transmitted  through  the  prism 
polarize  and  lost.)  A  final  lens,  L3,  focuses  the  light  onto  the  detector.  A 
pinhole  placed  just  in  front  of  the  detector  helps  reject  stray  light 
reflected  from  the  mount.  The  directional  coupler  is  simple  to  construct  and 
can  provide  uniform  excitation  of  all  bound  modes. 

Individual  bound  mode  groups  can  be  launched  in  a  step-index  waveguide 
using  a  combination  of  apertures  and  stops  in  the  spatial  filter  block  following 
lens  Lj .  Apertures  alter  the  size  of  the  beam  impinging  on  lens  L^.  This 
causes  the  launch  numerical  aperture,  NA^,  to  vary  according  to 

NAl  -  b/f  ,  (20) 
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Figure  4. 

Block  diagram  of  an  automated  OTDR  employing  a  Gian 
prism  polarizer  as  the  optical  directional  coupler. 

Only  the  horizontal  component  of  the  unpolarized,  fiber 
pigtail  output  is  launched  into  the  waveguide  under  test. 
Only  the  vertical  component  of  the  unpolarized  backscat- 
tered  signal  is  sent  to  the  detector.  Consequently,  the 
minimum  total  insertion  loss  of  this  coupler  is  6  dB. 
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where  2b  is  the  diameter  of  the  beam  entering  L2 ,  and  f  is  the  focal  length  of 
lens  L2 .  The  mode  groups  that  are  excited  can  be  determined  from 


m  _  1  b 

M  “  7TK  f  ’ 


(21) 


where  m  is  the  principal  mode  group  number,  and  M2  Is  the  number  of  bound  modes 
in  the  waveguide-5 

Although  polarization  rotated  by  90*  from  the  input  polarization  is 
detected,  this  usually  presents  no  measurement  problem.  In  standard  multimode 
glass  waveguides,  light  is  depolarized  within  a  few  meters.  This  distance  is 
small  compared  with  the  probe  pulse  width  of  the  OTDR;  consequently,  the  mea¬ 
surements  are  unaffected  by  a  polarization-sensitive  directional  coupler. 
Complications  occur  if  polarization-preserving  waveguides  are  to  be  measured. 
(Polarization-preserving  waveguides  include  single-mode  and  elliptical-core 
multimode  waveguides.)  Some  method  must  be  used  to  scramble  the  polarization 
if  all  modes  are  to  be  equally  and  independently  excited  in  these  waveguides. 

Operation  of  the  OTDR  is  straightforward.  System  synchronization  is 
provided  by  the  master  clock.  Timing  pulses  from  the  clock  provide  triggering 
information  for  both  the  transient  digitizer  and  the  pulse  circuitry  of  the 
laser.  The  detected  signal  is  amplified  and  fed  into  the  transient  digitizer. 

A  complete  OTDR  characterization  of  the  waveguide  Is  recorded  by  the  transient 
digitizer  for  each  laser  pulse.  Digitizer  resolution  Is  512  time  cells  by  512 
vertical  deflection  cells.  A  predetermined  number  of  individual  characteriza¬ 
tions  are  performed  and  averaged  to  Increase  the  signal-to-noise  ratio.  The 
resulting  waveform  is  sent  to  a  desktop  computer  for  processing  and  then  to 
the  X  Y  plotter  for  hard  copy  output. 

The  computer  program  which  oversees  data  acquisition  performs  three  func¬ 
tions.  It  determines  waveguide  length,  supervises  data  acquisition,  and 
analyzes  data  for  hard  copy  preparation.  Waveguide  length  is  determined  by 
measuring  the  elapsed  time  between  the  near-  and  far-end  Fresnel  reflections 
and  multiplying  by  the  group  velocity.  Data  reduction  consists  of  two  parts: 
calculating  the  average  optical  power  backscattered  from  each  point  along  the 
waveguide,  and  determining  the  local  attenuation.  The  result  is  a  pair  of  X  Y 
plots.  The  first  is  a  serailogarithmic  plot  of  backscattered  optical  power 
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versus  position  along  the  waveguide  (referred  to  as  the  power  plot).  The 
second  is  a  plot  of  local  slope  of  the  power  plot  (referred  to  as  the  loss 
plot).  The  loss  plot  is  obtained  from  the  raw  data  by  performing  a 
least-squares  fit  to  an  exponential  over  a  time  interval  equivalent  to  the 
probe  pulse  width.  The  least-squares  fit  is  performed  for  all  the  data  points 
in  the  data  array. 

If  the  waveguide  under  test  contains  substantial  attenuation  variations, 
the  computer  program  characterizes  the  waveguide  from  each  end  and  then  separ¬ 
ates  geometric  variation  from  actual  power  decay.  Local  waveguide  attenuation 
is  determined  from  a  least  squares  fit  to  the  data  of  the  actual  power  decay 
curve. 

Figure  5  is  a  characterization  of  a  waveguide  obtained  with  this  second 
program.  The  top  of  the  figure  shows  that  the  waveguide  exhibits  very  little 
geometric  variation  over  length.  Geometric  variation  is  found  from 

t(z)  -  10  log  .  (22) 

Three  separate  curves  are  shown  in  the  center  of  Figure  5.  The  solid  curves 
are  the  characterizations  from  end  1  and  end  2.  The  dashed  line  Is  the  actual 
power  decay  curve.  Notice  that  a  small  point  insertion  loss  (0.17  dB)  is 
present  400  ra  down  the  waveguide.  The  point  insertion  loss  is  even  more  pro¬ 
nounced  in  the  local  attenuation  plot  at  the  bottom  of  the  figure.  The  origin 
of  this  point  insertion  loss  is  a  single  microbend.  The  other  variations  in 
the  local  attenuation  are  diameter  variations. 

Key  performance  specifications  of  this  Instrument  are  listed  in  Table  1. 
The  directional  coupler  provides  a  mount  scatter  rejection  of  at  least  20  dB 
at  0.9  pm.  The  directional  coupler  exhibits  a  total  insertion  loss  of  just 
under  10  dB  which  includes  coupling  losses  into  the  fiber  pigtail  as  well  as 
losses  due  to  the  polarization  properties  of  the  prism  polarizer.  Typically, 
100  mW  of  source  power  is  coupled  into  the  waveguide  under  test.  The  launch 
beam  has  a  maximum  0.2  NA  and  a  spot  size  of  50  pm.  Typical  pulse  width  is 
130  nsec,  which  translates  into  a  pulse  26  m  long  inside  the  waveguide. 
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Table  l.  Key  Performance  Specifications  of  the  Automated  OTDR 


Parameter 

Typical  Performance 

Peak  power  coupled  into  waveguide 
at  0.9  pm 

100  mW 

Launch  NA 

0.2 

Launch  beam  spot  size 

50  pm 

Total  insertion  loss  of  direc¬ 
tional  coupler 

<10  dB 

Pulse  width 

130  nsec 

Spatial  resolution 

26  m 

Accuracy  of  length  measurements 

±2.5% 

Dynamic  range 

50  dB,  i.e.,  instrument  can  see 
through  a  waveguide  exhibiting 

25  dB  attenuation 

Fresnel  reflection  isolation  of 
directional  coupler 

Exceeds  20  dB 

Precision 

±0.023  dB/km  on  waveguide 
exhibiting  7.64  dB/km 
attenuation 

Accuracy 

Within  4% 

Length  measurements  are  accurate  to  about  2.5%,  limited  by  the  resolution  of 
the  transient  digitizer  and  the  accuracy  with  which  the  waveguide  group  veloc¬ 
ity  is  known. 

A  precision  test  consisting  of  25  successive  measurements  of  a  single 
waveguide  was  performed  using  this  apparatus.  The  waveguide  is  a  1-km-long, 
graded-index  waveguide  with  125  pm  o.d.,  62.5  pm  core  diameter,  and  0.18  NA. 
For  the  precision  test  the  probe  beam  spot  size  was  50  pm  with  a  launch  NA  of 
0.20.  The  mean  attenuation  was  7.64  dB/kra  at  0.9  pm  with  a  standard  deviation 
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of  ±0.023  dB/kra.  This  represents  the  highest  precision  yet  reported  for  an 
OTDR.  OTDR  accuracy  was  determined  by  comparing  OTDR  measurements  with  meas¬ 
urements  made  on  the  same  waveguide  using  the  conventional  cut-back  technique. 
A  total  of  16  cut-back  measurements  on  the  waveguide  yielded  a  mean  attenua¬ 
tion  of  7.81  dB/kra  with  a  standard  deviation  of  0.61  dB/km.  (This  measurement 
was  also  made  with  a  0.20  launch,  NA,  50-pm-dtameter  spot  size.)  The  agree¬ 
ment  between  measurements  is  within  0.17  dB/km  (<4%). 
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SECTION  4 


CONSIDERATIONS  OF  OTDR  SENSITIVITY  AND  REPEATABILITY 


As  we  pointed  out  earlier,  one  of  the  disadvantages  of  the 
backscattering  technique  is  that  the  exact  nature  of  the  backseat  ter ing  is 
unknown.  For  example,  to  obtain  accurate  attenuation  data,  it  is  necessary 
that  the  backscattered  power  be  strictly  proportional  to  the  power  present  In 
the  waveguide.  If  the  backscattered  power  is  a  constant  fraction  of  the  power 
in  the  waveguide  and  the  geometry  of  the  waveguide  does  not  change,  then  we 
have  the  simple  case  that  the  OTDR  signal  is  directly  proportional  to  the 
power  in  the  waveguide  (see  Equation  (10)).  It  was  first  assumed, &  and  later 
experimentally  proven,7  that  the  backscattered  power  has  a  A-4  dependence  and 
hence  is  Rayleigh  scattering.  Rayleigh  scattering  depends  only  on  the  mate¬ 
rial  composition  and  fabrication  technique  and  has  a  nearly  isotropic  angular 
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distribution  over  the  NA  within  the  waveguide.  Therefore,  in  Rayleigh- 
scattering-limited  waveguides,  the  backscattered  power  is  strictly  propor¬ 
tional  to  the  power  present  in  the  waveguide.  The  lowest-loss  waveguides 
manufactured  today  are  Rayleigh-scattering-liraited  waveguides.  Hence,  the 
OTDR  is  expected  to  perform  well  on  these  waveguides. 

An  additional  extremely  Important  question  is  the  distribution  of  the 
backscattered  power  among  the  bound  modes  of  the  waveguide. 

Knowledge  of  the  backscattered  distribution  is  important  because  it  has 
a  direct  impact  upon  OTDR  perturbation  sensitivity.  We  seek  to  locate  and 
characterize  perturbations  using  the  backscattering  method.  Since  each  mode 
propagates  with  a  characteristic  attenuation  (which  Is  modified  by  the  pres¬ 
ence  of  a  perturbation),  complete  characterization  of  the  perturbation 
requires  a  knowledge  of  both  the  launch  and  backscattered  modal  distributions. 
As  an  example,  consider  a  perturbed  simple  two-mode  waveguide  where  mode  1 
propagates  with  low  attenuation  and  mode  2  with  high  attenuation.  An  OTDR 
characterization  of  this  waveguide  could  result  in  four  different  permutations 
of  the  measured  attenuation,  depending  upon  modal  excitation. 
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Our  analysis  provides  an  equation  which  can  be  used  to  determine  modal 
excitation  of  backscattered  light  in  a-profile  waveguides.  We  assume  no  mode 
mixing  and  no  differential  modal  attenuation.  Under  these  conditions  a  WKB 
analysis9  indicates  that  the  radial  variation  of  power  is 


2,  \  2 

n  (r)  -  tu 

p(r)  =  — - -  ,  (23) 

no  "  "2 

if  all  modes  are  excited  equally  at  launch.  The  quantities  in  Equation  (23) 
are 
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(24) 


(25) 


The  local  numerical  aperture,  A(r),  is 


A(r)  *  n(r)  sin0c(r)  =*  [n2(r)  -  n2]1^2  ,  (26) 

where 

n  =  on-axis  refractive  index 
o 

n2  =  clad  refractive  Index 
a  =  core  radius 

a  =  a  constant 

Since  we  have  assumed  no  mode  mixing  and  no  differential  mode  attenuation,  the 
power  at  any  point  along  the  waveguide  is 
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where  PD  is  the  total  power  launched  Into  the  waveguide,  and  v  is  the 
attenuation  coefficient. 


The  fraction  of  light  which  is  scattered  in  an  annulus  of  radius  r  which 
is  trapped  in  backward  traveling  rays  is 
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where  we  assume  the  light  is  scattered  isotropically.  Consequently,  the 
radial  distribution  of  the  backscattered  power  is 
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where  S  is  the  fraction  of  the  pulse  which  is  scattered  at  z. 

The  mode  group  m,  excited  by  a  ray  which  enters  the  waveguide  at  radius 
r,  propagating  at  an  angle  0  with  respect  to  the  optic  axis,  can  be  found 
from5 
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where  M2  Is  the  number  of  modes  In  the  a-proflle  waveguide.  Inside  the 
w.i v»*gu i do  a  hound  ray  will  result  If  the  ray  forms  an  angle  so  that 


0  ‘  6  .  (31) 

v  n(r) 

Consequently,  the  mode  groups  which  are  excited  by  light  scattered  In  an 
annulus  of  radius  r  are 

»  (crt-2)/2a 

<r> :  4 

2  * 
n  (r) 


which  can  be  rewritten  as 


(32) 


Equation  (31)  provides  us  with  knowledge  of  which  mode  groups  are  excited  but 
does  not  determine  how  much  power  backseat tered  at  radius  r  Is  coupled  Into 
each  bound  mode.  To  determine  the  power  coupled  Into  each  bound  mode  we 
assume  that  each  mode  receives  equal  power  and  count  the  number  of  modes. 

Each  mode  group,  m,  Is  m-fold  degenerate.  The  lowest  order  mode  group 
which  Is  excited  by  light  scattered  from  an  annulus  of  radius  r  is  D-fold 
degenerate,  where 


,(u+2)/2 

°  =  M  (a ) 

The  next  higher  order  group  is  D  +  1  fold  degenerate;  the  highest  order  mode 
group  is  M-fold  degenerate. 
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Therefore,  the  total  number  of  excited  modes  is 


N  -  2  [D  +  (D  +  1)  +  ...  +  M]  «  M(M  +  1)  -  D(D  -  1)  ,  (34) 

where  the  factor  of  2  accounts  for  the  two  possible  polarizations.  Finally, 
the  power  back-scattered  from  radius  r  into  any  allowed  mode  is 


P  Sn2A  1  -  (r/a)a  2e  Y* 

raode^  *  2  M(M  +  1)  -  D(D  -  1)  *  U 

Figure  6  shows  the  relative  back-scattered  power  in  various  a-profile  wave¬ 
guides  as  a  function  of  mode  group  number.  Figure  7  shows  the  fractional 
change  in  mode  group  excitation  assuming  all  modes  are  equally  excited. 


Source  Considerations 


Once  the  data  acquisition  system  was  verified  we  began  looking  into  the 
question  of  the  optimum  source  for  this  application.  The  backscatterlng  which 
occurs  in  the  fiber  is  described  by  a  spatial  function,  b(z).  (We  assume  the 
fiber  is  infinitely  long  for  simplicity.)  If  a  small  intrusion  perturbation, 
f(z),  occurs  along  the  fiber,  the  total  backscattered  signal  is  written  b(z)  + 
f(z).  The  OTDR  probe  pulse  is  a  function  of  time,  t,  and  position,  p(ct  -  z). 
The  temporal  dependence  of  the  backscattered  signal,  B(t),  is 


B(t ) 


f- 


f 


(b(z)  +  f(z))  p(ct  -  z)  dz  =  (b(z)  p ( ct  -  z)  dz  +  lf(z)  p(ct  -  z)  dz 


/■ 


The  first  term  in  Equation  (36)  is  the  backscatterlng  of  the  unperturbed 
fiber.  The  second  terra  is  the  change  in  the  backscatterlng  caused  by  the 
intrusion  perturbation.  We  wish  to  maximize  the  second  term  by  tailoring  the 
shape  of  the  probe  pulse. 

Since  the  functions  in  Equation  (36)  are  convolutions,  the  sensitivity 
is  maximized  by  matched  filtering.  Ideally,  we  would  tailor  the  probe  pulse 
so  that  it  is  a  matched  filter  for  whatever  perturbation  is  present. 
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Figure  6. 

Relative  backscattered  power  in  a-profile  waveguides  as  a 
function  of  mode  group,  normalized  to  peak  power  and  back- 
scattered  by  a  step-index  waveguide. 
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Unfortunately,  there  are  a  large  number  of  potential  Intrusions  and  we  have  no 
prior  knowledge  of  which  intrusion  would  be  employed.  The  ability  to  generate 
a  matched  filter  for  each  one  would  cause  an  inordinate  level  of  complexity  in 
the  OTDR  source  electronics;  therefore,  a  simpler  approach  must  be  found. 

A  better  approach  for  a  real  system  is  to  employ  a  wide  probe  pulse. 
Without  resorting  to  detailed  mathematical  argument,  several  important  reasons 
for  using  a  wide  pulse  are  given.  If  the  probe  pulse  is  narrower  than  the 
perturbation,  the  sensitivity  is  less  than  optimum.  If  the  probe  pulse  is, 
roughly  speaking,  just  as  wide  as  the  perturbation,  the  sensitivity  is 
maximum;  however,  spatial  resolution  suffers.  For  near-term  applications  it 
is  imperative  to  know  that  an  Intrusion  is  present  rather  than  to  know  the 
fine  structure  in  the  backscattered  return.  We  conclude  that  to  maximize  the 
sensitivity  of  the  OTDR  for  a  variety  of  perturbations  a  wide  pulse  is  best. 

There  is,  therefore,  a  tradeoff  between  the  requirements  of  high  spatial 
resolution  and  high  OTDR  sensitivity.  To  increase  spatial  resolution  requires 
that  the  probe  pulse  width  be  decreased.  (Eventually,  in  the  quest  for  reso¬ 
lution  the  probe  pulse  width  is  reduced  below  the  width  of  the  perturbation 
and  sensitivity  suffers.)  Assuming  constant  peak  pulse  power,  a  shrinking 
pulse  width  reduces  the  backscattered  signal  power  proportionately. The 
required  detector  bandwidth  is  inversely  proportional  to  pulse  width;  hence, 
detector  noise  is  Inversely  proportional  to  pulse  width.  Signal-to-noise 
ratio  therefore  drops  as  the  square  of  the  pulse  width;  hence,  increasing 
resolution  reduces  sensitivity.  Sensitivity  may  be  increased  by  increased 
Integration  time.  However,  if  the  probe  pulse  width  is  reduced  by  a  factor  F, 
the  required  Integration  times  will  increase  by  F2. 

Laser  diode  sources  alone  are  incapable  of  the  peak  powers  required  for 
submeter  resolution.  Increasing  the  spatial  resolution  of  the  OTDR  requires 
either  an  Increase  In  the  peak  optical  power  in  each  pulse,  or  electrically 
removing  random  noise.  Since  we  are  already  pulsing  our  semiconductor  laser 
sources  at  their  maximum  forward  drive  current,  these  sources  require 
Increased  integration  times  to  Increase  spatial  resolution.  An  alternative  is 
to  go  to  other  laser  types  such  as  an  Nd:YAG  laser. 


Several  benefits  accrue  from  an  NdtYAG  source.  First,  there  is  more 
than  ample  peak  pulse  power.  In  fact,  most  of  the  available  power  would  be 
sacrificed  to  prevent  stimulated  effects  and  laser  damage  from  occurring  in 
the  waveguide.  In  addition,  the  1.06  and  1.3  pm  wavelengths  of  this  source 
are  much  closer  to  the  lowest  loss  window  in  current  fibers,  thus  permitting 
observation  of  substantially  greater  distances  with  the  OTDR. 

Several  problems  are  associated  with  this  solution.  As  the  pulse  width 
Is  reduced,  the  pulse  repetition  frequency,  prf,  goes  down.  A  reduced  prf 
requires  more  expensive  detection  and  signal  processing  electronics,  and 
requires  longer  data  acquisition  times.  There  is  also  a  question  of  the  reli¬ 
ability  of  the  source.  The  only  suitable  commercially  available  laser  which 
we  have  located  is  a  Quantel  laser.  One  of  these  lasers  is  already  owned  by 
HRL  and  is  notorious  for  its  alignment  and  reliability  problems.  This  laser 
requires  daily  realignment,  a  serious  drawback.  The  laser  also  exhibits  a 
short  mean-time  to  failure.  These  problems  suggest  that  the  best  initial 
choice  for  a  high  spatial  resolution  pTDR  is  a  semiconductor  laser  source  and 
integration  of  the  backscattered  signal. 
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SECTION  5 


OPTICAL  FIBER  PERTURBATION  MEASUREMENTS 

A.  DETECTION  OF  SPLICES  AND  PERIODIC  MICROBENDING 

Experiments  were  performed  to  determine  the  OTDR  signatures  of  splices 
and  periodic  microbending.  Waveguides  number  1  and  2  from  a  commercially 
available  Corning  10  fiber  cable  were  used  in  the  experiment.  Both  waveguides 
exhibit  essentially  constant  loss.  Fiber  1  exhibits  about  6.5  dB/km,  while 
fiber  2  exhibits  about  8.5  dB/km.  Figure  8  presents  the  OTDR  return  obtained 
after  fusing  the  two  waveguides.  The  Insertion  loss  of  the  splice  was  found 
to  be  0.2  dB  from  the  power  plot.  Notice  that  this  splice  was  obtained  by 
fusing  end  2B  to  end  IB.  The  loss  plot  shows  a  monotonic  decrease  in  the  loss 
in  fiber  1  following  the  splice.  The  monotonic  decrease  was  seen  in  each  of 
three  fusion  splices  made  with  these  waveguides  and  is  attributed  to  the  decay 
of  leaky  modes  generated  by  the  fusion  splice.  Finally,  notice  that  fiber  2 
now  exhibits  only  6  dB/km  average  attenuation,  compared  to  8.3  dB/km  before 
splicing.  This  gross  deviation  in  attenuation  is  not  yet  understood.  The 
only  known  change  that  occurred  in  the  experiment  is  that  the  launch  end  has 
been  exchanged.  We  are  now  launching  Into  end  2A,  whereas  the  Initial  charac¬ 
terization  was  performed  while  launching  into  end  2B.  A  variation  as  large  as 
2.3  dB/km  attenuation  when  launch  ends  ate  exchanged  has  not  been  observed 
previously. 

A  second  set  of  experiments  was  performed  to  attempt  to  isolate  a  signa¬ 
ture  for  periodic  raicrobending  structure.  For  this  experiment,  a  microbending 
loss  modulator  was  employed.  The  modulator  uses  a  series  of  rods  to  bend  the 
waveguide.  The  rods  are  spaced  3.5  mm  apart,  as  shown  in  Figure  9.  The  wave¬ 
guide  follows  a  serpentine  path  between  the  rods  exhibiting  a  bending  period 
of  7  mm.  Coupled  mode  theory  predicts  that  this  period  is  adequate  to  cause 
all  lower  order  modes  to  couple,  but  insufficient  to  cause  coupling  among  the 
highest  order  modes.  Mass  is  applied  to  the  upper  plate  in  50  gin  increments, 
yielding  the  following  peak-to-peak  amplitude  distortions  In  the  fiber  axis. 
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Figure  8.  OTDR  Characterization  of  a  fuse-spliced  fiber. 
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Figure  9.  A  microbeading  transducer. 
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Figure  10  is  a  typical  backscattered  signal  which  resulted  when  the  loss 
modulator  was  placed  immediately  after  the  splice.  Notice  the  exponential 
drop  in  attenuation  following  the  splice  and  modulator.  This  drop  is  due  to 
the  attenuation  of  leaky  modes  excited  at  the  splice.  Tt  was  observed  that 
the  loss  in  fiber  1,  immediately  adjacent  to  the  splice,  monotonically 
increases  with  distortion  amplitude;  howeve.  »  the  time  the  pulse  arrives  at 
the  far  end  of  the  waveguide,  the  loss  is  bacrt  i_  i  its  initial  (unperturbed) 

6.5  dB/kra  loss. 

The  modulator  was  also  placed  immediately  before  the  splice  and  the 
experiment  repeated.  The  results  were  nearly  identical.  Both  the  splice  and 
the  loss  modulator  introduce  a  point  insertion  loss  and  cause  mode  conversion 
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Figure  10.  0.9  pm  OTDR  characterization  of  a  fuse-spliced  fiber 

with  a  loss  modulator  placed  irameidated  after  the  splice. 
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to  leaky  modes.  If  the  two  perturbations  occur  in  series,  the  order  has 
little  effect. 

A.  final  experiment  was  performed  using  only  the  bending  loss  modulator 
and  fiber  2.  The  bending  loss  modulator  was  placed  immediately  after  the 
launch.  Figure  11  shows  the  result  of  this  experiment.  Note  that  the  change 
in  the  OTDR  characteristic  is  only  about  half  of  the  change  which  was  observed 
in  the  original  experiments  during  which  a  splice  was  present.  This 
difference  is  due  to  mode  conversion  in  the  splice. 

Extremely  small  displacements  of  a  microbend  transducer  have  been 
reliably  located  and  measured  with  the  OTDR.  The  transducer  consists  of  two 
mating  plates  with  milled  grooves  (1.5  mm  spacing)  into  which  stainless  steel 
rods  (500  pm  diameter)  are  cemented.  The  transducer  thus  provides  a  minimum 
mechanical  distortion  period  of  1.5  mm  which  can  be  Increased  by  removing  pins 
or  by  altering  the  angle  between  the  fiber  axis  and  the  pin  axes.  The  fiber 
is  subjected  to  a  total  of  thirty-one  bends  in  the  transducer. 

The  results  of  experiments  performed  on  this  transducer  are  shown  in 
Table  2.  Displacements  were  measured  using  a  linear  voltage  displacement 
transducer  calibrated  at  1.22  mV/ura.  The  small  displacements  measured  during 
the  experiments  with  the  hydrophone  transducer  lead  to  questions  of  accuracy, 
it  seems  unlikely  that  the  hydrophone  transducers  is  mechanically  perfect.  It 
is  particularly  unlikely  that  the  tops  of  the  pins  in  the  transducer  all  lie 
in  a  single  plane.  It  is  far  more  likely  that  if  a  plane  were  drawn  through 
the  mean  pin  height,  some  pins  would  protrude  above  that  plane  and  some  would 
be  recessed.  To  insure  that  the  tops  of  all  the  pins  lie  in  the  same  plane 
requires  some  sort  of  polishing  operation  be  performed  on  them. 

A  simple  alternative  approach  is  a  single  bend  transducer  employing  only 
three  pins.  If  a  thin  conductive  coating  is  also  applied  to  the  waveguide,  a 
conductivity  measurement  makes  it  possible  to  determine  when  this  fiber  is 
just  in  contact  with  all  three  pins.  This  is  a  convenient  reference  point 
corresponding  to  the  point  where  the  waveguide  just  begins  to  bend. 
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Figure  11.  Local  fiber  attenuation  immediately  following  the  splice  as  a 
function  of  the  peak-to-peak  amplitude  distortion  of  the  fiber 
axis. 
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Table  2.  Attenuation  of  the  Microbend  Transducer 


Experimental  Conditions 

Spat ial 
Period 

A,  mm 

Peak-to-Peak 

Distortion 

Amplitude, 

pm 

Two-Pass 
Insert  ion 
Loss,  dB 

Dry,  all  modes  excited 

1.5 

3 

0.30 

equally 

6 

0.62 

8 

1.16 

Dry,  all  modes  excited 

2.1 

10 

0.04 

equally 

13 

0.04 

14 

0.10 

20 

0.31 

Wet,  index  matching 

2.1 

0 

0.02 

fluid  n  =  1.51,  all 

5 

0.21 

modes  excited  equally 

11 

0.93 

B.  DETECTION  OF  A  SINGLE  LOOP  IN  A  FIBER 

Experiments  designed  to  collect  signatures  of  various  perturbations  have 
been  carried  out.  One  of  the  perturbations  Is  a  simple  fiber  loop. 

There  Is  a  simple  method  for  determining  the  presence  of  a  loop  using 
the  OTDR.  It  Is  known  that  modes  with  a  principal  mode  number,  ra, 
approximately  equal  to  or  greater  than  a  critical  value,  n^,  defined  by 


m 

c 


M 


» 


(37) 


are  essentially  and  Instantly  lost  in  a  simple  bend.  The  variables  In  this 
equation  are  the  core  radius,  a,  the  Index  difference  between  core  and  clad, 
A,  the  radius  of  the  bend,  R,  and  the  total  number  of  mode  groups,  M.  In 
step-index  waveguides  there  Is  a  simple  relat  "  ship  between  the  principal 
mode  number  and  launch  angle  of  a  collimated  beam.  That  relationship  is 


ra  _  sln8 
M  sinS 

max 


(38) 
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where  0  Is  the  angle  of  the  beam  with  respect  to  the  optic  axis,  and 
®max  Is  the  critical  angle  for  total  internal  reflection.  Hence,  we 
have  a  simple  method  of  determining  whether  a  fiber  has  a  bend  in  it,  and  even 
for  determining  what  the  radius  of  the  bend  is.  If  the  OTDR  measures  no 
excess  loss  at  the  location  of  the  loop  when  6=0,  and  large  excess  loss  when 
0  =  Omax,  a  loop  is  present.  (This  conclusion,  of  course,  assumes  that 
little  or  no  mode  mixing  occurs  in  the  fiber.)  If  we  measure  the  critical 
angle,  0C,  for  which  the  OTDR  begins  to  measure  large  excess  loss,  we  can 
deduce  the  bend  radius  from  the  expression 


A  standard  Corning  step-index  waveguide  was  used  for  this  experiment. 

The  waveguide  has  an  NA  of  0.184,  an  outer  diameter  of  125  urn,  and  a  core 
diameter  of  85  pm.  Most  of  the  waveguide  is  on  the  standard  12-cm  Corning 
shipping  spool.  Approximately  125  ra  have  been  wound  onto  an  identical  12-cm 
spool  (see  Figure  12).  A  loop  was  formed  between  the  two  spools  and  it  was 
placed  in  an  Index  matching  fluid.  The  initial  loop  diameter  was  12.6  cm. 
measurements  were  made  at  10  loop  diameters  between  2.1  and  12.6  cm.  The 
results  of  the  experiment  are  shown  in  Figure  13.  For  this  experiment  all 
mode  groups  were  excited  with  approximately  the  same  power. 

We  have  also  measured  the  excess  attenuation  induced  by  a  loop  using 
five  different  launch  modal  configurations.  The  configurations  are  summarized 
in  Table  3.  The  results  are  plotted  in  Figure  14,  along  with  a  simple 
theoretical  prediction  determined  from  the  formula: 


Figure  12.  The  experimental  setup  for  detecting  a  single  fiber  loop. 

A  single  loop  of  fiber  from  the  middle  of  a  long  fiber  is 
immersed  in  index  matching  fluid. 


LOOP  INSERTION  LOSS,  dB 
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Figure  13.  The  observed  attenuation  of  a  single  loop  of  waveguide 
when  all  modes  are  equally  excited. 
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where  R  Is  the  radius  of  the  loop.  The  curve  labeled  "A"  in  Figure  14  pres¬ 
ents  data  obtained  from  two  different  experiments.  For  the  first  experiment, 
the  focused  spot  was  placed  at  the  center  of  the  core.  For  the  second  experi¬ 
ment,  the  focused  spot  was  placed  at  the  core-clad  boundary.  Virtually  iden¬ 
tical  results  were  obtained  as  expected. 

The  larger  excess  loss  for  lower  order  mode  excitation  as  compared  to 
the  intermediate  mode  group  is  not  understood.  Note  that  the  OTDR  can  only 
detect  loops  as  large  as  8  cm  in  diameter.  Larger  diameter  loops  might  be 
located  if  the  signal~to-noise  ratio  is  improved. 

C.  MEASUREMENT  OF  THE  INSERTION  LOSS  OF  A  SINGLE  MICROBEND 

A  possible  technique  for  taping  an  optical  fiber  is  rale  robe nding. 
Microbends  are  microscopic  bends  at  the  optical  axis  which  are  short  enough  to 
cause  mode  conversion.  If  mode  convertlon  is  present,  power  which  was  bound 
to  the  core  can  be  converted  to  even  high  order  modes  until  it  is  finally  lost 
to  radiation.  The  power  which  is  lost  is  then  available  for  use  by  the 
intruder. 

Under  this  contract  we  have  measured  the  sensitivity  of  the  OTDR  to  the 
presence  of  the  microbends.  Single  miccobends  as  small  as  10  vm  in  amplitude 
and  1.5  mm  in  extent  were  observable.  The  OTDR  response  was  found  to  be 
nearly  quadratic  in  raicrobend  amplitude. 

The  computerized  OTDR,  described  in  Section  3,  was  used  to  characterize 
the  microbend. 


Table  3.  Details  of  the  Launch  Modal  Configurations  used 

in  the  Initial  Loop  Detection  Experiments 


Launch 

Configuration 

Aperture 

Opening, 

mm 

Stop 

Diameter, 

mm 

Launch 

NA 

Focused 
Spot  Size, 
pm 

Approximate 
Mode  Groups, 
Launched 

A 

10.0 

0 

0.20 

5.0 

0-53 

B 

6.5 

0 

0.13 

7.6 

0-34 

C 

3.5 

0 

0.07 

14.0 

0  -  18 

D 

6.5 

3.5 

0.07  <  NA  <  0.13 

— 

18  -  34 

E 

10.0 

6.5 

0.13  <  NA  «:  0.20 

— 

34  -  35 
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Figure  14.  F.xcess  loss  as  a  function  of  inverse  loop  radius  for  a  step 
Index  fiber. 
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A  300-m-long  high-quality  step-index  waveguide  was  used  Cor  this 
experiment.  The  core  diameter  is  84  pm,  the  outer  diameter  is  126  pm,  the 
numerical  aperture  is  0.20,  and  the  mode  mixing  length10  is  670  m. 

The  raicrobending  experimental  setup  is  shown  in  Figure  15.  Sixty  meters 
of  the  fiber  are  wound  on  spool  A  and  the  remainig  270  m  are  wound  on  spool  B. 
Between  the  spools  is  the  raicrobend  transducer,  which  is  in  turn  mounted  on  a 
precision  vertical  translation  stage  (Klinger  UH-80PP).  This  stage  provides  1 
±  0.5  pm  vertical  steps  on  command  from  the  controller.  The  transducer  plates 
have  many  parallel  grooves,  spaced  1.5  mm  apart,  milled  into  the  contacting 
surfaces.  Two  steel  drill  blanks,  0.53  ram  in  diameter,  are  cemented  into  the 
upper  plate  in  adjacent  slots.  One  drill  blank  is  cemented  into  the  bottom 
plate  in  a  slot  centered  between  the  two  upper  rods.  This  arrangement 
provides  a  single  microbend  of  1.5  mm  period  when  the  waveguide  axis  is  normal 
to  the  drill  blank  axis.  (The  period  can  be  changed  by  rotating  the 
transducer).  An  integrating  sphere  collects  the  0TDR  probe-pulse  power 
transmitted  by  the  waveguide.  A  lockin  amplifier  (PARC  124)  measures  the 


Figure  15.  The  raiccobending  experimental  setup. 


power  collected  within  the  sphere.  Simultaneous  measurement  of  the 
transmitted  power  provides  an  Insertion-loss  measurement  check  on  the  accuracy 
of  the  OTDR  data. 

Preliminary  experiments  demonstrated  that  the  OTDR  is  extremely  sensitive 
to  the  presence  of  ralcrobending  perturbations.  This  forced  the  use  of  a  posi¬ 
tive  method  for  determining  when  the  pins  just  touch  the  waveguide  The  solu¬ 
tion  chosen  was  to  immersion  silver  the  waveguide.  Approximately  1  pm  of 
silver  is  deposited  on  the  waveguide  by  this  process.  The  overcoated  wave¬ 
guide  is  placed  in  the  microbending  fixture  and  electrical  leads  are  attached 
to  all  three  drill  blanks.  The  bottom  plate  of  the  microbending  fixture  is 
stepped  up  in  1  pm  steps  to  the  point  where  electrical  continuity  is  just 
obtained  between  all  three  pins.  Microbend  amplitude  is  measured  with  respect 
to  this  reference  point. 

Typical  OTDR  microbend  characteristics  are  presented  in  Figure  16.  The 
top  half  of  the  figure  presents  characterizations  of  the  unperturbed  waveguide 
and  microbends  of  10-  and  25-pm  amplitude.  The  bottom  half  of  the  figure 
presents  the  change  in  the  local  attenuation  attributable  to  the  presence  of 
the  ralcrobend.  These  difference  curves  were  obtained  by  subtracting  the  local 
slope  of  the  unperturbed  waveguide  characterization  from  the  local  slope  of 
the  characterizations  obtained  when  the  three  bending  pins  just  touched  the 
waveguide,  when  a  10-pm  perturbation  was  present,  and  when  a  25-pm  perturba¬ 
tion  was  present.  The  insertion  loss  of  the  microbend  was  determined  by  inte¬ 
grating  the  area  under  the  difference  curves.  Note  that  there  is  no  reflec¬ 
tion  occurring  at  the  microbend.  Only  loss  is  present. 

The  experimental  results  are  shown  in  Figures  17  and  18.  Figure  17  shows 
insertion  loss  as  a  function  of  microbend  amplitude  for  two  different  rod 
spacings,  1.5  mm  and  2.1  mm.  The  solid  curves  were  obtained  from  the  OTDR 
data,  while  the  dotted  curves  are  insertion  loss  measurements.  The  maximum 
difference  between  corresponding  measured  values  is  0.074  dB.  The  hard  pri¬ 
mary  coating  was  removed  from  a  small  length  of  the  fiber.  Adjacent  segments 
of  coated  and  uncoated  waveguide  were  immersion  silvered.  Microbends  were 
Introduced  Into  each  segment  and  Identical  insertion  losses  were  measured. 

The  data  fit  well  to  a  square  law  formulation,  as  shown  in  Table  4.  (The 
single  data  point  at  25-pm  distortion  amplitude,  2.1-mm  period  was  deleted 
when  calculating  the  best  power  law  fit.) 
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Figure  16.  Typical  OTDR  microbend  characteristics. 
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Figure  17.  Microbead  Insertion  loss  dependence  on 
the  square  of  the  ralcrobend  amplitude. 
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A  =  1 .5  mm 

COATED  WAVEGUIDE 


A  dependence  on  the  square  of  the  distortion  amplitude  is  expected.  The 
average  power  coupling  which  occurs  between  modes  of  any  principal  mode  group 
is  proportional  to  the  power  spectrum  of  the  fiber  distortion.^  The  distor¬ 
tion  applied  to  the  fibers  is  a  truncated  sinusoid  with  amplitude  A.  Thus  the 
bending  power  spectrum  is  proportional  to  the  square  of  the  distortion 
amp l itude. 

Only  a  slight  dependence  on  launch  mode  configuration  was  observed  (see 
Figure  18).  This  figure  compares  the  measured  results  for  the  cases  of  low, 
low  and  intermediate,  and  all  modes  launched.  The  waveguide  has  a  maximum 
mode  group  number  of  48  at  the  measurement  wavelength  of  0.9  pm.  The  low  mode 
launch  excites  mode  groups  1  through  20..  The  low  and  intermediate  mode 
launch  excites  groups  1  through  38. 

Relative  insensitivity  to  mode  group  excitation  is  expected  since  this  is 
a  step-index  waveguide.  In  step-index  waveguides,  the  difference  in  propaga¬ 
tion  constant  between  adjacent  modes  depends  upon  the  mode  group  number  M 
according  to1'’ 


where  A  is  the  relative  index  difference  between  core  and  clad,  a  is  the  core 
radius,  and  M2  is  the  number  of  bound  modes  in  the  waveguide.  Coupling  will 
occur  between  mode  groups  i  and  j  if  their  propagation  constants  are  related 
to  the  mechanical  period.  A,  of  the  perturbation  through  the  expression 


(42) 


Perturbations  with  spectral  components  of  longer  or  shorter  period  will 
not  couple  power  efficiently.  The  largest  intraraode  group  jump  occurs  between 
the  highest  and  second  highest  mode  groups.  No  significant  conversion  of 
hound  mode  power  to  radiation  modes  can  occur  unless  power  is  coupled  across 
this  gap.  Nonsinusoidal  mechanical  perturbations  short  enough  to  couple  these 
mode  groups  can  be  expected  to  contain  sufficient  low-frequency  spectral  com¬ 
ponents  to  couple  all  lower  lying  mode  groups.  Longer  period  perturbations 
tmpart  insufficient  momentum  to  the  mode  group  to  permit  coupling  across  the 
largest  "propagation  constant  gap." 
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Table  4.  Experimental  Results 


Induced  Attenuation,  dB 

A  =  1 . 5  mm 

A  =  2.1  mra 

Distortion 

0TDR 

Insertion  Loss 

OTDR 

Insertion  Loss 

Amplitude, 

pm 

Data 

Data 

Data 

Data 

0 

0 

0 

0 

0 

5 

— 

— 

- 

— 

10 

0.026 

0.026 

— 

— 

15 

0.058 

0.058 

— 

— 

20 

0.091 

0.104 

— 

— 

25 

0.155 

0.176 

0.079 

0.005 

30 

0.232 

0.271 

0.083 

0.037 

35 

0.343 

0.388 

0.115 

0.054 

40 

0.470 

0.524 

0.160 

0.086 

45 

0.648 

0.693 

— 

— 

50 

0.901 

0.893 

- 

— 

Best  fit 

power 

law 

2.2 

2.2 

2.3 

2.9 

Coefficient 
of  correla¬ 
tion 

0.99 

1.00 

1.00 

0.99 
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D.  MEASUREMENTS  OF  FIBER  DIAMETER  VARIATIONS 

If  all  the  high  order  modes  of  a  multimode  waveguide  are  populated  with 
an  equal  amount  of  power,  and  if  the  diameter  of  the  waveguide  decreases,  then 
some  of  these  modes  are  cut  off  and  the  power  they  carry  is  lost.  Virtually 
all  waveguides  have  some  diameter  fluctuations.  Power  is  available  to  an 
intruder  if  it  can  locate  a  diameter  variation  point  and  if  power  is  lost  at 
this  point. 

The  OTDR  can  be  used  to  Identify  points  where  such  an  intrusion  might  be 
attempted.  Consider  a  narrow  pulse  of  light  propagating  along  a  waveguide. 

The  pulse  suffers  a  position  dependent  attenuation  described  by  the  attenua¬ 
tion  coefficient,  Yf(z),  composed  of  a  scattering  term,  Yfs(z),  and 
an  absorption  term,  Yfa(z).  The  change  in  pulse  power,  dPf(z),  is 
linearly  proportional  to  the  pulse  power  according  to 

dPf(z)  =■  ~Yf(z)  Pf(z)  dz  ,  (43) 


which  can  be  rewritten  as 


Pf(z)  =  exp 


z 

J  Yf(x)  dx 
o 


(44) 


where  PQ  is  the  Initial  peak  pulse  power.  The  backscattered  power,  P^,  is 
some  backscattering  function,  f(z),  of  the  power  at  position  z.  The  backscat- 
tering  function  depends  upon  waveguide  geometric  properties  (NA  and  core  diam¬ 
eter)  as  well  as  the  scattering  coefficient.  Scattering  sites  encountered  by 
the  probe  pulse  are  most  commonly  microscopic  Rayleigh  scattering  sites. 

These  sites  scatter  nearly  isotropically,  permitting  a  simple  first  order 
calculation  by  the  waveguide  acceptance,  NA,  at  the  scattering  site  by  4x. 

The  backscattering  function  is  simply  the  fraction  multiplied  by  the  scatter¬ 
ing  coefficient.  For  a  step-index  waveguide  the  backscattering  function 
is10 


f  (z) 


2 


NAjLzJ. 

n2(z) 


(45) 
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The  backscattered  signal  suffers  position  dependent  attenuation,  Yb(z),  on 
the  return  trip.  Collecting  these  results,  we  find  that  the  power  which  <s 
backscattered  from  the  pulse  at  point  z  and  is  subsequently  detected  is 


P(z)  =  CP^fCz)  exp 


|“  J  (Yf(x)  +  Yb(x))dx 


The  quantity  C  is  the  coupling  efficiency  of  the  directional  coupler.  Now 
consider  a  waveguide  containing  core  diameter  fluctuation. 

Since  the  number  of  bound  modes,  N,  depends  on  the  core  radius,  a, 
according  to 


fka  NA) 


changes  in  radius  cause  corresponding  changes  in  the  number  of  bound  modes. 
Mode  coupling  cannot  occur  if  the  period.  A,  is  large.  Higher  order  modes 
must  become  leaky,  causing  excess  attenuation  which  can  be  lumped  into  the 
absorption  coefficient.  If  the  bound  mode  excitation  is  constant,  Independent 
of  position  and  independent  of  direction  of  travel,  then  a  sinusoidal 
variation  in  core  radius  of  the  form 


a(z)  =  aQ  [l  +  B  sin(£z)] 


modulates  the  attenuation  coefficient  according  to 


Y  -  Yt  +  Ya"i(z)  sin(£z)  , 


where  £  Is  2x/A,  and  m(z)  depends  upon  the  properties  of  the  waveguide.  If 
the  backscatter  function  and  the  modulation  depth  are  constant,  the  power 
backscattered  as  a  function  of  position  is  given  by 

P(z)  »  P  fe  2^Zexp  (-  (1  -  cos(£z))|  .  (50) 


Figure  19  Is  the  first  reported  observation  of  diameter  variations  made  with 
the  OTDR.11  The  figure  shows  a  characterization  of  a  step-index  waveguide 
which  was  manufactured  without  the  aid  of  a  diameter  stabilization  system. 

The  lack  of  stabilization  resulted  in  peak  diameter  variations  of  ±32  pm 
around  the  mean  waveguide  diameter  of  120  pm.  The  period  of  the  variations  is 
about  100  m.  Note  the  one-to-one  correspondence  between  peaks  in  backseat- 
tered  signal  and  local  minima  in  the  waveguide  diameter. 

Much  more  subtle  diameter  variations  can  be  seen  in  most  commercially 
available  waveguides.  The  signature  of  one  such  waveguide,  a  high-quality, 
graded-index  waveguide  manufactured  by  ITT,  is  shown  in  Figure  20.  The  signa¬ 
ture  shown  here  was  obtained  when  all  modes  were  equally  excited  by  the  probe 
pulse.  There  are  many  very  small  deviations  from  linearity  in  this  character¬ 
ization.  To  accentuate  them,  the  local  slope  of  the  curve  was  measured  and  is 
also  shown  here.  Notice  that  the  local  attenuation  varies  nearly  sinusoidally 
with  position. 

It  was  not  possible  to  remove  the  primary  coating  from  this  fiber  in 
order  to  measure  the  actual  outer  diameter.  (This  particular  waveguide  was 
destined  for  an  application  requiring  that  its  initial  high  strength  be  main¬ 
tained.  The  fiber  had  passed  a  200  kpsi  proof  test.)  To  verify  that  diameter 
fluctuations  were  responsible  for  this  observed  attenuation  behavior  we  con¬ 
tacted  the  manufacturer  and  asked  for  their  diameter  chart.  ITT  graciously 
provided  that  chart.  The  measured  fiber  diameter  is  shown  in  Figure  21.  Note 
that  there  are  sinusoidal  diameter  variations;  however,  the  amplitude  is  quite 
small,  about  2.5  pm  peak  corresponding  to  an  RMS  variation  <1.5%  of  the  0D. 
Also  shown  is  the  local  attenuation  from  Figure  20  and  the  superposition  of 
the  two  curves.  There  Is  excellent  correlation  between  them.  From  this  data 
we  estimate  the  sensitivity  of  the  OTDR  to  sinusoidal  diameter  variations  to 
be  about  1.5  dB/km/pm  change  in  core  diameter.  Note  also  the  high  spatial 
frequency  diameter  variations.  These  variations  are  also  visible  in  the  local 
waveguide  attenuation  data  In  spite  of  the  fact  that  these  variations  are 
typically  10  m  in  width  (which  is  much  narrower  than  the  OTDR  probe  pulse 
width  of  25  m). 
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Figure  19.  First  reported  observation  of  diameter  fluctuations  made 
with  the  0TDR. 
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To  try  to  improve  our  understanding  of  these  observations,  a  special 
uncoated  waveguide  was  drawn  (Figure  22,  curve  D.)  The  waveguide  has  a  nomi¬ 
nal  outer  diameter  of  126  pm.  This  diameter  was  held  constant  to  better  than 
±1/4  pm  for  the  first  150  m.  Then  sinusoidal  diameter  fluctuations  of  150  m 
period  and  10  pm  amplitude  were  drawn  into  the  fiber.  The  silica-germanium- 
phosphorus  core  of  this  st  p-index  fiber  has  a  diameter  of  70  pm.  The  silica- 
boron  buffer  layer  is  3  pm  wide.  The  clad  is  126  pm  in  diameter  and  is  made 
from  an  Hereaus-Amersil  T08-WG  tube.  The  NA  is  0.20. 

The  forward  OTDR  signature  is  shown  in  Figure  22,  curve  Pf.  For  this 
characterization,  a  50-pm  spot,  0.2  NA  launch  configuration  was  employed. 

Thus  all  mode  groups  were  excited  with  nearly  equal  power.  The  signature  is 
similar  to  the  diameter  variation,  but  there  are  unexpected  complications. 
First,  note  that  the  local  peaks  in  backscattered  signal  occur  when  the  fiber 
diameter  is  minimum  as  was  noted  previously.  However,  there  is  a  shift 
between  the  location  of  the  diameter  maxima  and  the  local  backscattering 
minima.  This  shift  is  equivalent  to  a  60°  phase  shift.  Note  also  that  the 
signature  looks  like  the  sum  of  a  linear  function  and  a  sinusoid.  Observa¬ 
tion:  if  linear  and  sinusoidal  functions  are  added,  the  peaks  of  the  compos¬ 

ite  shift  in  phase  with  respect  to  the  original  sinusoid  by  tsin**1  of  the 
slope  of  the  linear  function.  Also  note  that  the  average  slope  of  this  signa¬ 
ture  is  30°. 

Observe  that  whenever  the  fiber  diameter  decreases,  the  slope  of  the 
backscattered  return  also  decreases  and  can  become  positive.  Normally  this 
would  be  interpreted  as  a  decrease  in  the  local  attenuation.  However,  if  the 
modes  of  the  fiber  are  filled  at  the  launch,  and  if  the  diameter  is  reduced, 
then  the  number  of  bound  modes  are  converted  to  radiation  modes  and  loss  is 
incurred.  The  conclusion  to  be  reached  is  that  OTDR  characterizations  of 
waveguides  with  large  diameter  variations  must  be  interpreted  carefully. 

Tiie  waveguide  was  also  characterized  from  the  opposite  end  (Figure  22, 
curve  Pjj)*  Although  the  launch  conditions  are  the  same,  the  backscattered 
signal  is  different.  No  longer  is  there  a  backscattered  maxima  when  the  fiber 
diameter  Is  a  minimum,  however,  there  Is  a  slight  bump  on  the  signature  which 
does  correspond  to  diameter  minima.  Also  apparent  is  the  fact  that  the  local 
backscattering  minima  are  in  phase  between  the  two  characterizations.  The 
local  backscattering  maxima  are  out  of  phase.  The  only  significant  difference 
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between  these  two  characterizations  Is  that  In  the  first  case  the  diameter 
initially  decreased,  while  in  this  case  It  initially  Increases. 

DiVita  and  Rossi12  developed  a  method  for  separating  the  geometric  varia¬ 
tion  from  the  power  decay  (Figure  23).  They  observed  that  the  backscattered 
signal  consists  of  two  terms:  a  geometric  terra  and  an  exponential  term.  The 
exponential  term  contains  information  about  the  total  attenuation  coefficient, 
u.  The  geometric  term  is  proportional  to  the  forward  scattering  and  the  fiber 
cross-sectional  area.  That  is,  it  is  proportional  to  the  geometry  of  the 
fiber.  If  the  waveguide  is  characterized  from  the  forward  direction,  the 
resulting  backscattered  signal  is  given  by 

F  =  G(z)  exp(az)  .  (51) 

If  it  is  characterized  from  the  backward  direction,  it  is  given  by 

B  =  G(L  -  z)  exp  (a(L  -  Z)}  .  (52) 


If  the  geometric  terms  are  equal,  the  geometric  mean  of  the  forward  and  back¬ 
ward  characterizations  is  the  geometric  term,  and  the  square  root  of  the  ratio 
is  the  exponential  terra.  Thus, 


and 


G(z)  »  /FB  *  constant  x  G(z)  , 

P(z)  *  /F/B  51  constant  x  exp(az)  .  (53) 


Figure  22,  curve  G  shows  the  geometric  terra,  along  with  the  fiber  diam¬ 
eter.  The  geometric  terra  is  a  reasonably  good  copy  of  the  diameter  variation 
except  for  a  60°  phase  shift. 

Figure  22,  curve  P  shows  the  exponential  decay  term.  Only  a  small  devia¬ 
tion  From  linearity  is  apparent.  Taking  the  derivative  of  this  curve  provides 
the  local  attenuation  data  shown  in  Figure  22,  curve  L.  Again,  there  is  good 
correlation.  The  local  attenuation  is  maximum  when  the  diameter  is  minimum  as 
expected;  however,  the  60°  phase  shift  is  once  again  evident.  The  small 
amplitude  high  frequency  oscillations  are  caused  by  quantization  error  in  the 
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transient  digitizer.  This  same  experiment  was  repeated  for  two  other  launch 
modal  configurations  with  identical  results.  These  launch  configurations  were 
filling  the  lowest  1/3  of  the  mode  volume  and  the  highest  1/3. 

An  independent  measurement  of  the  scattered  power  was  also  made.  The 
experimental  setup  Is  shown  in  Figure  24.  A  He-Ne  laser  beam  was  chopped  and 
focused  onto  the  core  of  the  waveguide.  The  launch  NA  was  0.2.  A  cladding 
mode  stripper  was  employed  immediately  after  the  launch.  Light  scattered  at 
right  angles  to  the  fiber  axis  was  detected  by  a  PIN  photodiode  after  passing 
thcough  a  2-mm  aperture.  A  lock-in  amplifier  measured  the  detected  signal. 

The  waveguide  was  wound  one  layer  deep  on  a  30-cm  diameter  drum.  Position 
along  the  waveguide  was  determined  by  multiplying  the  meter  stick  reading  by 
the  winding  pitch.  The  experimental  results  are  shown  in  Figure  25.  The 
average  attenuation  of  the  waveguide  was  14  dB/km  at  0.63  pm.  Each  data  point 
was  corrected  to  account  for  this  mean  attenuation  before  it  was  plotted.  The 
scattered  power  is  180°  out  of  phase  with  the  diameter  variations  as  expected. 
Power  was  launched  into  both  ends  with  qualitatively  identical  results. 

In  conclusion,  diameter  variations  can  be  detected  with  the  OTDR.  Fea¬ 
tures,  in  general,  are  not  exactly  anticorrelated,  but  exhibit  some  phase 
shift.  The  backscattered  signal  is  complex,  but  is  apparently  relatively 
independent  of  the  launch  modal  configuration.  The  sensitivity  of  the  OTDR  to 
diameter  variation  Is  about  1.5  dB/km/ pm  change  in  core  diameter.  Diameter 
changes  as  small  as  1  pm  have  been  detected. 

E.  PRECISION  PERIODIC  MICROBEND  EXPERIMENT 

Microbend  transducer  is  a  potential  method  for  tapping  power  out  of  opti¬ 
cal  fibers.  The  mechanism  for  bending-induced  loss  in  an  optical  fiber  is  the 
transfer  of  guided-raode  power  to  the  radiation  modes.  A  distortion  of  spatial 
periodicity  small  enough  to  cause  mode  conversion  is  called  a  microbend. 
Furthermore,  a  periodic  disturbance  of  periodicity,  ft,  may  couple  two  modes 
labeled  l  and  j  If 


2* 

ft 


> 
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Figure  24.  Experimental  setup  for  scattering  measurement. 
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Figure  25.  Measured  scattering  power  along  the  fiber  with  diameter 
variation. 
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where  0^  and  0j  are  the  longitudinal  propagation  constants  of  modes  i  and 
j.  Determining  those  values  of  £2  that  can  be  expected  to  cause  mode  conver¬ 
sion  is  then  reduced  to  a  calculation  of  63 i j • 

For  a  fiber  with  the  refractive-index  profile  given  by 


2,  \  2 

n  (r)  =  nQ 


r  is  the  radial  coordinate,  nQ  =  n(0)  is  the  Index  at  core  center,  a  is  the 
core  radius,  A  is  the  fractional  difference  in  refractive  index  between 
cladding  and  core  center,  and  a  is  a  constant  having  the  value  a  =  2  for  a 
parabolic  index  fiber  and  a  =  "  for  step-index  fibers.  The  longitudinal 
propagation  constants  are  given  by 
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m 


n  k 
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2ot/(a+2) 


1/2 


(56) 


where  k  =  2n/X  (X  is  the  optical  wavelength  in  vacuum),  m  is  the  mode  label, 
and  M,  the  maximum  value  of  m,  is 


M  =  an  k 
o 


(57) 


The  last  two  equations  are  valid  for  highly  multimode  fibers,  and  in  this 
limit  the  total  number  of  modes  is  equal  to  M2  . 

Defining  60m  =  0^  -  0m, 


,  /  A  \1 /2  /  \a-2/a+2 
60  ( —  j 

m  a  \ a  +  2/  \M/ 


(58) 


For  a  parabolic  index  fiber  (a  =  2),  the  last  equation  yields  the  interesting 
result  that  mode  group  spacing  is  independent  of  mode  number: 


60 


(59) 


The  spatial  wavelength  of  Interest  for  mode-(..>upltng  considerations  is 

thus 
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(60) 


o  2jr_  2na 

C  "  66  "  (2M1/2  * 

For  the  Corning  1150  graded-index  fiber  (2a  =  65  pm  and  A  ■  0.008),  the  mode 
coupling  period  is  =*  1.6  mm. 

For  step-index  fibers  (a  =  °°),  60m  becomes 


The  largest  value  of  60  occurs  for  m  »  M,  the  smallest  value  of  60  for  m  *  1: 


max  a 


.  _  2A1/2  1  _  2 

rain  a  M  ,  2 
n  ka 
o 


For  the  Corning  1028  step-index  fiber  (2a  ”  85  pm,  A  *  0.008,  iIq  •  1.46), 
the  corresponding  spatial  wavelengths  are 


min  60 


- =  cm  . 

max  60 

min 
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From  these  numerical  examples  we  note  that  distortions  of  ft  ~  10  cm  will  mix 
low-lying  modes,  but  smaller  ft's  are  required  to  mix  the  higher  order  modes.  For 
the  parabolic-index  fiber,  simple  bending  couples  only  nearest-neighbor  mode 
groups.  Thus  only  a  single  ft  is  required  to  cause  significant  loss.  On  the 
other  hand,  since  no  simple  selection  rule  exists  for  the  step-index  fiber, 
all  wavelengths,  ft  <  ®min>  are  expected  to  produce  significant  loss.  No 
single  wavelength  will  cascade  power  from  low-lying  modes  of  the  step  fiber  to 
the  radiation  modes  as  effectively  as  does  ftc  for  the  parabolic  fiber. 

The  preceding  analysis  leads  to  an  interesting  conclusion:  if  a  pure 
sinusoidal  perturbation,  with  spatial  period,  ft  <  ftmjn,  is  generated  in 
a  step-index  fiber,  then  it  is  possible  to  create  mode  mixing  of  two  specific 
mode  groups  without  net  power  loss. 

In  order  to  verify  this  point  we  carried  out  the  following  experiment. 

The  setup  is  almost  Identical  to  that  used  to  measure  the  Insertion  loss  of  a 
single  microbend.  The  only  difference  is  the  microbend  transducer  is  composed 
of  two  corrugated  plates  with  precise  groove  spaclngs. 

A  photomask  was  designed  and  procured.  The  pattern  contains  a  periodic 
array  or  dark  strips  100  pm  in  width  with  1.2  mm  period.  Two  silicon  wafers 
(2  in.  diameter)  with  polished  (100)  surfaces  were  sputtered  with  a  Si02 
film.  This,  in  turn,  is  used  as  the  mask  for  the  subsequent  preferential 
chemical  etching  of  silicon.  A  sketch  of  the  profile  of  the  finished  silicon 
wafers  is  shown  in  Figure  26.  These  wafers  are  mounted  on  quartz  flats, 
which,  in  turn,  are  mounted  on  vacuum  knuckles  that  allow  the  two  plates  to  be 
in  perfect  parallelism  when  pressed  against  each  other. 

The  top  plate  is  stationary  while  the  bottom  plate  is  on  a  vertical 
translation  stage  driven  by  a  digitally  controlled  stepping  motor  which  has  a 
resolution  of  ~1  pm. 

A  spool  of  multimode  step-index  fiber  ~260  meters  long  was  used  in 
this  experiment.  The  microbend  transducer  is  located  about  half  way  along  the 
fiber.  Figures  27(a)  and  (b)  shows  an  0TDR  plot  of  the  fiber  subject  to 
microbending  of  various  fiber  perturbation  amplitudes.  The  fiber  sees  about 
40  periods  of  microbend.  It  is  seen  that  a  displacement  amplitude  as  small  as 
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SILICON  SUBSTRATE 


1  Pm  is  detectable  on  the  differential  attenuation  .lot.  In  this  case  all 
modes  are  launched.  We  have  also  selectively  launched  only  lower  order  modes 
and  only  higher  order  modes.  The  OTDR  plots  of  these  three  different  cases 
did  not  exhibit  significant  differences  as  one  might  expect.  We  suspect  the 
fiber  might  not  be  mode-preserving  since  it  was  wound  on  a  small  drum  without 
special  precautions.  Unfortunately,  the  silicon  plates  were  damaged  during 
one  of  the  runs,  precluding  further  investigations. 


DIFFERENTIAL  ATTENUATION.  dB/km 


SECTION  6 


CONCLUSIONS 


The  OTDR  has  been  shown  to  be  a  unique  instrument  capable  of 
nondestructively  performing  rapid,  accurate,  precise  measurement  of  most  opti¬ 
cal  waveguide  parameters.  It  is  the  only  instrument  that  provides  information 
on  the  length-dependence  of  waveguide  properties.  Total  attenuation,  scatter¬ 
ing  loss,  absorption  loss,  location  and  insertion  loss  of  defects,  insertion 
loss  of  connectors,  splice  loss,  mode-mixing  length,  and  the  effect  of 
environmental  perturbations  can  all  be  measured.  If  a  waveguide  lias  substan¬ 
tial  defects,  characterization  of  the  waveguide  from  both  ends  allows  the 
spatial  variation  of  waveguide  geometry  to  be  separated  from  the  actual  power 
decay. 

A  series  of  experiments  have  been  performed  which  place  the  OTDR  on  a 
firm  experimental  foundation.  Experiments  demonstrate  that  great  precision 
can  be  achieved  in  backscatter  measurements.  The  assumption  of  Rayleigh  scat¬ 
tering  has  been  verified.  The  tremendous  power  of  the  OTDR  to  locate  and  char¬ 
acterize  perturbations  has  been  demonstrated.  Diameter  fluctuations,  micro- 
bends,  splices  and  macrobend  (e.g.,  loop)  are  among  the  perturbations  which 
have  been  measured  with  the  instrument. 

Theoretical  considerations  indicated  that  maximum  perturbation  sensitiv¬ 
ity  of  the  OTDR  is  achieved  by  Increasing  the  source  power  until  the  threshold 
for  stimulated  effects  is  reached.  This  maximizes  sensitivity  regardless  of 
the  relative  shapes  of  the  probe  pulse  and  the  perturbation  response.  Tailor¬ 
ing  the  probe  pulse  shape  to  be  a  matched  filter  for  the  perturbation  results 
in  a  further  increase  in  sensitivity  at  the  cost  of  greatly  increased  complex¬ 
ity  In  the  source. 

Several  characteristics  of  perturbation  signatures  have  been  found. 
Perturbations  which  are  much  smaller  in  spatial  extent  than  the  probe  pulse 
are  characterized  by  the  shape  of  the  probe  pulse.  If  the  probe  pulse  width 
is  much  smaller  than  the  spatial  extent  of  the  perturbation,  the  OTDR  response 
is  the  impulse  response  of  the  perturbation.  Fresnel  reflection  sites  can  he 
separated  from  the  Rayleigh  scattered  signal  by  reducing  the  probe  pulse 
width.  Reducing  the  probe  pulse  width  decreases  the  Rayleigh 


scattered  signal.  Fresnel  reflection  depends  only  upon  an  Index  discontinuity 
across  a  dielectric  interface;  hence,  reducing  the  probe  pulse  width  does  not 
affect  the  amplitude  of  the  Fresnel  reflection. 

Experiments  on  splices  indicate  that  it  is  very  difficult  to  splice 
together  two  waveguides  so  perfectly  that  the  OTDR  is  unable  to  detect  the 
splice.  In  every  observed  case,  an  insertion  loss  was  measured  and  conversion 
of  bound  mode  power  to  leaky  modes  was  observed. 

Experiments  on  measuring  the  excess  loss  of  a  single  fiber  loop  indicated 
that  the  OTDR  technique  is  capable  of  detecting  loops  as  large  as  8-cm  in. 
diameter.  Larger  diameter  loops  can  be  located  if  the  system  signal-to-noise 
ratio  is  improved. 

Large  as  well  as  small  amplitude  microbends  in  a  step-index  waveguide  can 
easily  be  located  and  measured  with  the  OTDR.  The  microbend  transducer  intro¬ 
duces  an  insertion  loss  and  causes  mode  conversion  to  leaky  modes.  Under  this 
program  we  have  detected  the  presence  of  a  single  microbend  less  than  10  pm  in 
amplitude  and  the  presence  of  a  periodic  microbend  (~40  periods)  of  amplitude 
less  than  2  pm. 

Further  research  is  necessary  in  order  to  understand  the  modal  dependence 
of  the  OTDR  perturbation  signatures. 
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